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ABSTRACT 

The angular-diameter distance Da of a galaxy cluster can be measured by combining its 
X-ray emission with the cosmic microwave background fluctuation due to the Sunyaev- 
Zeldovich effect. The application of this distance indicator usually assumes that the cluster is 
spherically symmetric, the gas is distributed according to the isothermal /3-model, and the X- 
ray temperature is an unbiased measure of the electron temperature. We test these assumptions 
with galaxy clusters extracted from an extended set of cosmological N-body/hydrodynamical 
simulations of a ACDM concordance cosmology, which include the effect of radiative cool- 
ing, star formation and energy feedback from supernovae. We find that, due to the tempera- 
ture gradients which are present in the central regions of simulated clusters, the assumption of 
isothermal gas leads to a significant underestimate of Da- This bias is efficiently corrected by 
using the polytropic version of the /^-model to account for the presence of temperature gra- 
dients. In this case, once irregular clusters are removed, the correct value of Da is recovered 
with a 5 per cent accuracy on average, with a ~ 20 per cent intrinsic scatter due to cluster 
asphericity. This result is valid when using either the electron temperature or a spectroscopic- 
like temperature. When using instead the emission-weighted definition for the temperature 
of the simulated clusters. Da is biased low by ^ 20 per cent. We discuss the implications 
of our results for an accurate determination of the Hubble constant Hq and of the density 
parameter We find that, at least in the case of ideal (i.e. noiseless) X-ray and SZ obser- 
vations extended out to t^qq, Hq can be potentially recovered with exquisite precision, while 
the resulting estimate of ^Im, which is unbiased, has typical errors Ailm — 0.05. 

Key words: large-scale structure of Universe - galaxies: clusters: general - cosmology: mis- 
cellaneous - methods: numerical 



1 INTRODUCTION 

The Sunyaev-Zeldovich (SZ) effect JSunvaev & Zeldovich|[T97l 
is the distortion of the Cosmic Microwave Background (CMB) 
spectmm due to the scattering of the CMB photons off a popula- 
tion of electrons. At radio frequencies, the typical size of this dis- 
tortion for a thermal distribution of electrons with temperature of 
about 10 keV is at the level of 10~*. This effect has now been de- 
tected for a fairly large number of clusters of gal axies ("e.g. lRephaelil 
Il995t lBirkinshawlll999l : ICarlstrom et"n]l200l for reviews). For 
almost three decades it has been recognized that the combination 
of X-ray and SZ observations of galaxy clusters provides a direct 
measurement of the cosmic distance scale, under the assump tion of 
spherical symmetry for the intra-cluster gas distribution (e.g. | Guii d 
ll978tlSilk & WhitJl978tlCavaUere et alJl979HBirkinshawll979h 



The method is based on the different dependence on the electron 
number density, rie, of the X-ray emissivity (oc n^Tl^'^ for ther- 
mal bremsstrahlung; here Te is the electron temperature) and of the 
SZ signal (oc UeTe). 

Due to the crucial role played by the assumption of spher- 
ical symmetry, a great deal of efforts have been spent either to 
select individual clusters hav ing ve ry relaxed and regula r mor- 
phology (e.g . Holzaofel et Zlll997l: iH ughes & Birkinshawj ri998t 
lOrainge et alJl2002l : iBonamente et al. 2004), or to build suitable 
cluster samples over which avera ging out the uncertainti es due 
to intrinsic cluster ellipticity (e.g. [Mason e tal||2001 : Ree se et alj 
l2002HUdomDrasert et aljl2004ljones et alJl2005il) . These analyses 
have provided estimates of the Hubble constant, Hq, which are gen- 
erally consistent with those obtained from the Cepheid distance 
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scale (e.g. iFreedman et alj|200l|) or inferred from the spectrum 
of the CMB anisotropics (e.g. Isoergel et alj|2003h . although with 
fairly large uncertainties. Although the dominant source of uncer- 
tainty is probably represented by the contaminat ion of the SZ s ignal 
by the CMB and point-sources (e.g. Udompras ert et all2004) . sig- 
nificant errors are also associated to cluster asphericity, clumpy gas 
distribution and incorrect modeling of the thermal structure of the 
intra-cluster medium (ICM). 

So far, the limited number of high-redshift clusters with both 
SZ and X-ray observations, with their relatively large uncertain- 
ties, made the calibration of the cosmic distance scale mostly sen- 
sitive to the value of Hq, while no significant constraints have 
been placed on the values of the matter density parameter Q.m 
and the c osmological const ant. In the co ming years, ongoing X - 
ray (e.g.. lMullis etaPlZOOSh . optical (e.g.. lGladders & Yeell2005h . 
and planned or just started SZ surveys^ promise to largely increase 
the number of observed clusters out to 2 ~ 1.5. This may well open 
the possibility to use SZ/X-ray cluster observations to place con- 
straints on the Dark Matter and Dark Energy content of the universe 
jMolnar et all2002h . This highlights the paramount importance of 
having observational uncertainties under control. 

In this respect, numerical hydrodynamical simulations of 
galaxy clusters may offer an important test-bed where to quan- 
tify observational biases and keep the corresponding uncertainties 
under control. For instance, eliminating from the SZ and X- 
ray signal leaves a sensitive dependence of the angular-diameter 
distance, Da, on the electron temperature (see §2). On the other 
hand, temperature measurements of the ICM have been so far en- 
tirely based on fitting the X-ray spectrum to a suitable plasma 
model. How close is the resulting spectral temperature to the elec- 
tron temperature depends on the complexity of the thermal struc- 
ture of the ICM (e.g.. ,.Mazzotta et al. 2004). Hydrodynamical sim- 
ulations of clusters offer a natural way to quantify the bias intro- 
duced by replacing the electron temperature with the X-ray tem- 
perature. Furthermore, the standard assumption in the SZ/X-ray 
calibration of the cosmic distance scale is that of the isothermal 
ICM, while X-ray observations of clusters clear ly show the pres- 
ence of significant temperature gradients (e.g. [ M arkevitc h et alj 
ll998tlDe Grandi & Molendl20oilvikhlinin et all2005h . To o~ 
come this problem, several authors estimate the bias introduced 
by the isothermal approximation, finding t hat the distance can 
be biased by ^ 2 per cent (e.g. ■Bir kirishaw' & Hughes! 1 1 9941 : 
lUdomprasert et al J [2004 : Hol zapfel et alj Il997h . Simulations of 
galaxy clusters naturally produce temperature gradients that, at 
least at large r adii, are clos e to the o bserved ones (e.g., Loken et al. 
l2002HBorgani et al. 2004: IKav et ab 2004). Therefore, simulations 
can be used to quantify the bias introduced by the assumption of 
isothermal gas. Finally, using a representative set of galaxy clusters 
in a cosmological framework allows one to trace the distribution 

^ See, for example, the dedicated interferometer arrays: 

• AMI: http : / / www . mrao . ca^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 

• AMiBA: httpj^^www^^as^^^^^^^^^^^^^^^^^^^^ 

• SZA: h ttp: //astro ■ uchicago .edu/szel 

or the bolometers: 

• ACBAR: |http : / / cosmology .berkeley ■ edu/group/swlh/ acbar/| 

• ACT jhttp : / /www .h ep . upenn . edu/ ^a ngelica/ act /act .html| 

• APEX:^^^^^^^^^^^berkeley^^^du^a^^^^ 

• 01iinpQ: |http : //oberon . romal ■ infn. itTl 

• Planck: |http : / /astro ■ estec.esa.nl/Planck/| 

• SPT jhttp :/ /astro . uchicago . edu/ spt7| 



of ellipticity and, therefore, to calibrate the corresponding scatter 
in the measurement of the distance scale. Nowadays, cosmologi- 
cal hydrodynamical codes have reached a high enough efficiency, 
in terms of both achievable resolution and description of the gas 
physics, to provide a realistic description of the processes of for- 
mation and evolutio n of galaxy clusters (e.g. 'Borgani et alj |2004t 
iKravtsov et all2003) . For instance, Kazantzidis et al. ( 2004) found 
that halos in hydrodynamical simulations including cooling are sig- 
nificantly more spherical than in non-radiative simulations. Since 
the assumption of sphericity is at the basis of the X-ray/SZ method 
to estimate Da, this highlights the relevance of properly model- 
ing the physics of the ICM for a precise calibration of the cosmic 
distance scale. 

The purpose of this paper is to understand the impact of the 
above discussed systematics on the calibrations of the cosmic dis- 
tance scale from the combination of SZ and X-ray observations, 
by analyzing an extended set of hydrodynamical simulations of 
galaxy clusters. These simulations have been performed using the 
TREE-l-SPH GADGET-2 code I Sp ringel et al-2001. : ..Springel200j) . 
for a concordance ACDM model, and include the processes of ra- 
diative cooling, star formation and supernova feedback. The set 
of simulated clusters contains more than 100 objects having virial 
masses in the range (2-20) X lO^^h'^MQ. 

The plan of the paper is as follows. In Section 2 we review the 
method to estimate the angular-diameter distance from X-ray and 
SZ cluster observations in the case of a polytropic equation of state 
for the ICM, and discuss the different definitions of temperature 
that are used in the analysis of the simulations. In Section 3 we 
describe the set of simulated clusters and the procedure to generate 
X-ray and SZ maps. We present our results in Section 4, where we 
show the results on the accuracy of the measurement of Da- We 
discuss and summarize our main results in Section 5. 



2 Da from combined X-RAY AND SZ 
OBSERVATIONS 

The combination of the ICM X-ray emission with the SZ flux 
decrement provides a direct measure of the angular-diameter dis- 
tance. D a, of galaxy clusters (Silk & White 1978; Cavaliere et ^ 
Il97gt lBirkinshawlil9'79l) . The method takes advantage of the dif- 
ferent dependences of these two quantities on the electron number 
density. Tie (the former is oc while the latter is oc rtg). When 
combined, these two quantities provide the physical dimension of 
the cluster along the line-of-sight and, from the cluster angular size. 
Da, if the cluster is spherically symmetric. 

The Comptonization parameter y, as measurable from obser- 
vations of the SZ effect, is 



y 



ksTe 



(1) 



where fcs is the Boltzmann constant, gt is the Thompson cross 
section, is the mass of the electron, c is the speed of light and 
the integration is along the line of sight. By definition, it provides 
a redshift-independent measure of the total thermal content of the 
cluster. 

The X-ray surface brightness Sx is 



Sx = 



47r(l + z)4 



nenHKH(T)M 



(2) 



where uh is the hydrogen number density of the ICM, A.^h (T) is 
the cooling function (normalized to n^UH)- To compute the X-ray 
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emissivity of the simulated clusters, we assume the cooling func- 
tion taken from a Raymond-Smith code (Raymond & Smith 197% 
for a gas of primordial composition (Xh = 0.76 and Xhe ~ 0.24) 
for a fully ionized ICM. We compute the X-ray emissivity in the 
[0.5-2] keV energy band, which is used in several combined SZ/X- 
ray analyses relying on t he ROS AT-PSPC data for the X-ray imag- 
ing part (e.g.t Reese e t al. 2002; Jones et al. 2005). Using bands ex- 
tending to higher energies, as appropriate for Chandra and XMM- 
Newton observations, would produce no change in the final results. 



2.1 The polytropic /3-model 

A common procedure adopted to extract Da from the combination 
of eqs. Q and J2j is based on modeling the electron density profile 
with the /3-model, 



;(r) = HeO 



1 + 



-3,3/2 



(3) 



iCavaliere & Fusco-Femianol Il97d) . where Ueo is the electron 
number density in the cluster center, r is the distance from the clus- 
ter center, rc is the core radius and /3 is the power-law index. 

As for the temperature structure of the ICM, a number of anal- 
yses of X-ray data independently show that galaxy clusters are 
far from being isothermal. Significant negative gradients character- 
ize the temperature profiles of galaxy clusters, at least on scales 
R>0.2R2oo (e.g. Markevitch et al. 1998; De Grandi & Molendi 
[2002; Pr att & Arnau d 2002; P iffaretti et alJl2005l; IVikhlinin et all 
12005 ), with positive gra dients assoc iated only to the innermost 
cooling regions (e.g. All en et alj200ll) . The dynamic range covered 
by the SZ signal extends on scales which are relatively larger than 
those sampled by the X-ray emission. For this reason, one may ex- 
pect that a systematic effect is introduced by assuming the ICM to 
be isothermal when combining X-ray and SZ observations. Since 
the SZ signal has a stronger dependence on the ICM temperature 
than the X-ray one, the effect of assuming an isothermal ICM, in 
a regime where the temperature is decreasing, may lead to predict 
i/(S)-profiles which are shallower than the intrinsic ones. 

In order to account for the presence of temperature gradients, 
we introduce a polytropic equation of state, p oc p'' , which relates 
the gas pressure p to the density p, where 7 is the polytropic in- 
dex (7 = 1 for isothermal gas). The three-dimensional temperature 
profile is thus 



T4r) = T,_ 







V~' T 

1 = JeO 









-1 -3/3(7-l)/2 



(4) 



where Teo is the temperature at the cluster center. Using the above 
expression for the temperature profile in the definition of the Comp- 
tonization parameter of eq.Q gives 



yW = yo 



1+ TT 



2-1 (l-3/37)/2 



where the Comptonization parameter at the cluster center is 

fcsTeo ^r(3/37/2 - 1/2) 



yo = DArieoOcfTT- 



m,c- ■ r(3/37/2) 
Similarly, we obtain the X-ray surface brightness profile 

^ ^ ^ ^ {l-6/3[(7+3)/4]}/2 

Sx{d) ~ Sxo 

where the central surface brightness is 



1+^ 



(5) 



(6) 



(7) 



Sxo = DAn^QOc 



4^/7^ p.H 



r(3/3[(7 + 3)/4]-l/2) 
AMTeo) r(3/3[(7 + 3)/4] '^^^ 



In deriving the above equation, the dependence of the cooling func- 
tion on Te is assumed to be a power law, T", with index a = 0.5. 
This is valid in the case of pure bremsstrahlung emission and rep- 
resents a good approximation in the case of bolometric emissivity. 
However, our emissivity maps are build in the [0.5-2] keV band. 
In this energy range the cooling function is significantly flatter, due 
to the contribution of metal emission lines, which is relevant in the 
case of relatively cool systems {Te < 2 keV) (e.g. lEttoril200(]l) . In 
order to test the effect of approximating the cooling function with 
a bremsstrahlung shape, we repeated our analysis also in the bolo- 
metric band and found variations in the final distance estimates by 
< 10%. 

Finally, by eliminating rieo from eqs.j5| and Jsj, we obtain the 
angular-diameter distance 



Da 



A. 

Sxo 



ksTei 



AeH(T'eo)/ie/^H 1 
47r3/2cr|,(l + 2)'l 6»e 



r(3/37/2) 
r(3/3/2 - 1/2) 



r(3/3[(7 + 3)/4]-l/2) 
r(3/3[(7 + 3)/4]) 



(9) 



For 7 = 1, the above expression reduces to that usually adopted 
in observational analyses based on combin ing X-ray and SZ clus- 
ter observations (e.^. . Reese et al., 2002: Udomprasert et aLll2004l 
iBonamente et alj2004b . which relies on the assumption of isother- 
mal gas. 

It is worth reminding here that, while simulations are rather 
successful in reproducing th e negative te mperature gradient in 
the outer cluster regions (e.g. lEvrard 61^11996: Eke et alJll99^ 
iLoken et al J 120021 ; iRasia et al l2004i) ~thev generally produce too 
steep profiles in the central cluster regions, especially when radia- 
tive cooling is included (e.g. Katz & White 1993; Tornatore et sl\ 
i2003i;IValdamini 2003; Bo rgani et alj i2004). Since observed clus- 
ters are characterized by a core region which is closer to isother- 
mality than the simulated ones, we expect that the effect of using 
a polytropic temperature profile when analyzing simulated clusters 
is larger than the actual effect taking place in real clusters. 



2.2 Definitions of temperature 

In the observational determinations of Da from X-ray and SZ 
observations with eq.j9j one relies on the X-ray temperature ob- 
tained from the spectral fitting. Such a spectral temperature is gen- 
erally different from both the actual electron temperature, which 
appears in the expression of Da, and the emission-weighted tem- 
perature, which is often used as a proxy to the spectral tempera- 
turein the analysis of hydrodynamical simulations of clusters (e.g. 
lEvrardetal-1996') . 

If rie.i and Ti are defined as the electron number density and 
the temperature carried by the i~th simulation gas particle, then the 
electron temperature is defined by 



Te — ^ ^ f^e,i Ti / ^ ^ f^e.i , 



(10) 



which coincides with the mass-weighted temperature in the limit 
of a fully ionized plasma of uniform metallicity. Analogously, the 
emission-weighted temperature is 



Tew — 



E.Am; 



(11) 
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where the coohng function A(r) can be computed over an energy 
band, comparable to that where the X-ray spectrum is fitted in ob- 
servational data analyses. In the following, we compute the emis- 
sivity in the [0.5-7] keV band. 

However, Mathiesen & Evrard (2001) were the first to show 
that the emission-weighted temperature does not necessarily rep- 
resent an accu rate approximation to the spectroscopic temperature. 
[Mazzotta et al. (2004) have further motivated and quantified this 
difference, connecting it to a thermally complex structure of the 
ICM. These authors suggested an approximate expression for the 
spectroscopic temperature, the spectroscopic-like temperature: 

T.i = !' , (12) 



,a+3/2 



where a is a fitting parameter. iMazzotta et alJ i2004) have shown 
that eQ. <12> with a — 0.75 closely reproduces the spectroscopic 
temperature of clusters at least as hot as 3 keV, with a few per cent 
accuracy, after excluding all the gas parti cles colder than 0.5 keV 
from the sums in eQ. <12> . More recently, ( Vikhlinin'200^ has gen- 
eralized the above expression for T^i to include the cases of lower 
temperatures and arbitrary ICM metallicity. 

In the following, besides using the electron temperature, we 
also perform our analysis by relying on the temperature proxies 
of egs. jl It and <12> . Therefore, comparing the results based on 
the electron temperature and on the spectroscopic-like tempera- 
ture provides a check of the bias introduced by using the X-ray 
temperature in the estimate of Da, a bias possibly present also 
in the analysis of real data. Furthermore, the comparison between 
emission-weighted and spectroscopic-like temperature provides a 
hint on the bias introduced in the simulation analysis when using 
an inaccurate proxy to the X-ray temperature. It is worth reminding 
here that, due to the finite time for electron-ion thermalization, the 
corresponding electron and ion temperature may differ, for instanc e 
as a consequence of continuous shocks (e.g.. lYoshida et all2005l) . 
A sizable difference among these two temperatures may induce a 
bias in the estimate of the distance scale. 

Except for using different definitions of temperature, we do 
not investigate the effect of a realistic observational setup for the 
detection of both the SZ and X-ray signals. Besides the statisti- 
cal errors associated to time-limited exposures, we also neglect 
the effects of systematics (e.g., instrumental noise, foreground and 
background contribution from contaminants, etc.). A detailed anal- 
ysis of the contam inatio ns in the SZ signal has been provided by 
iKnox et all i2004) and bv lAghanim et alJi20o3h . A comprehensive 
description of the instrumental effects on the recovery of X-ray 
observables, calibrated on hydrodynamical si mulations, has b een 
provided by 'Gardini et al. (2004) (see also Rasia et al. 2006). In 
this sense, our analysis will be based on ideal maps, which are free 
of any noise. We defer to a future analysis the inclusion of the errors 
associated to realistic X-ray and SZ observational setups. 



3 THE SIMULATED CLUSTERS 

The sample of simulated galaxy clusters used in this paper has been 
extracted from the large-scale cosmological hydro-N-body simula- 
tion of a "concordance" ACDM model with Jim = 0.3 for the 
matter density parameter at present time, JIa = 0.7 for the cos- 
mological constant term, Jib ~ 0.019 h^^ for the baryons density 
parameter, h = 0.7 for the Hubble constant in units of 100 km 
s^^Mpc^^ and — 0.8 for the r.m.s. density perturbation within 
a top-hat sphere having comoving radius of 8 /i~^Mpc. We refer 



to lBorgani et alj i20o3) (B04 hereafter) for a detailed presentation 
of this simulation, while we give here only a short description. 

The run, performed with the massively parallel Tree-l-SPH 
code GADGET-2 (Sorineel et al. 2001; Soringel 2005), follows the 
evolution of 480'^ dark matter particles and an equal number of gas 
particles in a periodic cube of size 192/i~^ Mpc. The mass of the 
gas particles is rrigaa = 6.9 x 1O*^~^M0, while the Plummer- 
equivalent force softening is 7.5h~^ kpc at jz = 0. Besides gravity 
and hydrodynamics, the simulation includes the treatment of ra- 
diative cooling, the effect of a uniform time-dependent UV back- 
ground, a sub-resolution model for star formation from a multi- 
phase interstella r medium, as well as gala ctic winds powered by 
SN explosions ( Soringel & Hernauisll '2003). At z = we extract 
a set of 117 clusters, whose mass, as computed from a friends-of- 
friends algorithm with linking length fe = 0.15 (in units of the mean 
interparticle distance) is larger than lO^** h~^MQ. 

Due to the finite box size, the largest cluster found in the cos- 
mological simulation has Tc ~ 5 keV. In order to extend our anal- 
ysis to more massive and hotter systems, which are mostly rel- 
evant for current SZ observations, we include four more galaxy 
clusters having i\fvir > 10^^ h~^MQ^ and belonging to a differ- 
ent set of hydro-N-body simulations (Boraani et al. 2006). Since 
these objects have been obtained by re-simulating, at high resolu- 
tion, a patch of a pre-existing cosmological simulation, they have 
a better mass resolution, with mgas = 1.69 x W^h~^ Mq. These 
simulations have been performed by using the same code with the 
same choice of the parameters defining star-formation and feed- 
back. The cosmological parameters also are the same, except for a 
larger power spectrum normalization, erg = 0.9. 

Therefore, our total sample comprises 121 objects, spanning 
the range of spectroscopic temperatures Tsi ~ 1 — 9 keV, out of 
which 25 have Tsi > 2.5 keV and only four have Tsi > 5 keV. 
The corresponding temperature distribution is reported in FigureQ 
Quite apparently, our set of clusters on average samples a lower 
temperature range with respect to that covered by current SZ ob- 
servations. For this reason, we will discuss in the following the sta- 
bility of our results when selecting only the high end of the temper- 
ature distribution. Since our set of simulated clusters covers a rela- 
tively low temperature range, we can safely ignore any relativistic 
corrections to the SZ signal (e.g., Itoh et al. 1998). 

Around each cluste r we extract a s pheric al region extending 
out to 6 Rvir- Following IPiaferio et alji2005h . we create maps of 
the relevant quantities along three orthogonal directions, extending 
out to about 2 Ruir from the cluster center. Each map is a regular 
256 X 256 grid. 

In the Tree-l-SPH code, each gas particle has a smoothing 
length hi and the thermodynamical quantities it carries are dis- 
tributed within the sphere of radius hi according to the compact 
kernel: 



8 



M/(x) = — 2(1 



x> 1 



(13) 



where x — r/hi and r is the distance from the particle position. 
We therefore distribute the quantity of each particle on the grid 

Here and in the following, the virial radius, Rvir, is defined as the ra- 
dius of a sphere centered on the local minimum of the potential, containing 
an average density, Pvir, equal to that predicted by the spherical collapse 
model. For the cosmology assumed in our simulations it is pvir — lOOpc, 
being pc the cosmic critical density. Accordingly, the virial mass, Mvir, is 
defined as the total mass contained within this sphere. 
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Figure 1. The distribution of spectroscopic-like temperatures for tlie set of 
simulated clusters. The dashed and the solid lines are for the whole sample 
and for the subset of clusters classified as regular (see text), respectively. 
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Figure 2. Maps of the X-ray brightness and SZ y parameter for a regular 
simulated cluster having virial mass M^ir ~ 1.4 ■ lO^*/i~^M0, R500 = 
0.53h^^ Alps and Tgi = 2.2 keV. The map extends out to 2 i?viri so that 



it covers a physical scale of 6.05 h' 
done with a 256 X 256 pixelization. 



'^Mpc for this cluster Each map is 



points within the circle of radius hi centered on the particle. Specif- 
ically, we compute a generic quantity qjk on the {j, k) grid point 
as qjkdp = J q{r)dldp = ^ qi{mi/ pi)wi where dp is the pixel 
area, the sum runs over all the particles, and Wi oc J W{x)dl is the 
weight proportional to the fraction of the particle proper volume 
rrii/pi which contributes to the {j, k) grid point. For each particle, 
the weights are normalized to satisfy the relation '^Wk ~ 1 
where the sum is now over the grid points within the particle circle. 
When hi is so small that the circle contains no grid point, the parti- 
cle quantity is fully assigned to the closest grid point. Figures|2|and 
|3|show an example of the X-ray surface brightness and SZ maps 
and of the temperature maps of a relaxed cluster in our simulation 




Figure 3. Maps of the emission-weighted (top left), mass-weighted (top 
right) and spectroscopic-like (bottom left) temperature for the same cluster 
of Figurel2l 



(with M„,r ~ 1.4 • lO"/i"^M0), which we use in the following 
as an example. 



4 RESULTS 

In this section we present our results on the reliability of the usual 
procedure to recover the angular-diameter distance from the com- 
bination of the SZ and X-ray emission of clusters, by using both 
the isothermal and a more general polytropic equation of state for 
the ICM. Since the procedure to determine Da (z) is known to be 
particularly sensitive to the presence of cluster substructures and 
irregularities, the first step of our analysis is to select a subset of 
relaxed and regular clusters. Our selection criterion is based on vi- 
sual inspection of the X-ray and SZ maps, as well as on the profiles 
of the X-ray surface brightness and the Comptonization parameter. 
Although this is clearly not an objective criterion, it is quite similar 
to the criteria used to classify real clusters. Besides the example of 
a relaxed cluster shown in Figure |2| in Figure |4] we also show an 
example of a cluster that we classify as unrelaxed and, as such, is 
excluded from our analysis. Overall, we select 71 relaxed clusters 
from the initial sample of 121 simulated clusters (the four hottest 
clusters are all included in this subset). The temperature distribu- 
tion of these objects is represented in figureQwith a solid line. We 
have 19 clusters in this subset with Tsi > 2.5 keV. We extend our 
statistical sample by realizing "observations" of each cluster along 
three orthogonal lines of sight and treating them as three indepen- 
dent objects. 



4.1 Results from the isothermal model 

Unlike X-ray observations, current observations of the SZ effect 
in clusters do not allow to perform any spatially resolved analysis. 
For this reason, the commonly adopted procedure is to determine 
the parameters 9c and /3, which determined the /3-model density 
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Figure 4. The same as Fig|2] but for an unrelaxed simulated cluster having 



profile, from the X-ray imaging alone, along with the normaliza- 
tion Sxo- The SZ profile is then used to obtain the central value of 
the Comptonization parameter, yo, by using the /3-model parame- 
ters as determined from the X-ray profile. By following this proce- 
dure, we fitted all the profiles out to i?5oo , which is defined as the 
radius encompassing an average density of 500 times the critical 
cosmic density. We point out that iisoo corresponds to the typical 
outermost radius where X-ray observations provide surface bright- 
ness and temperature profiles. We exclude from the analysis the 
central regions of the clusters, within 1/20-Rvir, which are strongly 
affected by gas cooling and are close to the numerical resolution of 
the simulations. 

In Figure|5|we show the profiles of the X-ray surface bright- 
ness and Compton-j/ for the example cluster of Fig|2| along with 
the best-fitting /3-model for the isothermal case. For this relaxed 
cluster, the /3-model provides a rather good fit to the X-ray pro- 
file. Only the central point, which is anyway excluded from the fit, 
is higher than the /3-model extrapolation, as a consequence of the 
high-density gas residing in the cluster cooling region. The result- 
ing values of the fitting parameters for this cluster are /? = 0.835 
and rc/Rsoo = 0.196. Fitting the y{R) profile with eq.j5}, after 
setting P = 0.835 and 7=1, leads to an underestimate of j/o- 
In fact, the resulting Compton-j/ profile is significantly shallower 
than measured (Figure|5}- This result follows from neglecting the 
presence of negative temperature gradients. Consequently we tend 
to underestimate Da, because Da oc j/q [see eq.j9j]. In this par- 
ticular case we underestimate Da by 56 per cent. 

This result for one particular cluster is confirmed by the distri- 
bution shown in the left panel of Figure|6|(see also Table0. In this 
figure we report the distribution of the ratios D'^X'^ / D^a^'^ between 
the recovered and the true values of the angular-size distance. Such 
results clearly demonstrate that the angular-size distance is biased 
low, on average, by more than a factor two, as a consequence of the 
underestimate of yo induced by the assumption of an isothermal 
ICM. In order to verify a possible temperature dependence of the 
Da distribution, in the left panel of Fig|5|we also show the results 
for the clusters with Tsi in the range 2.5-5 keV and for those hot- 



ter than 5 keV. While the latter are too few to allow any meaningful 
conclusion, the clusters at intermediate temperature have a distribu- 
tion which is statistically consistent with that of the whole sample. 
This indicates the absence of any obvious trend of our results with 
the cluster size. The results reported in this figure have been ob- 
tained by using the electron-weighted temperature estimate for the 
simulated clusters. If we had used the emission-weighted temper- 
ature, we would have obtained an even stronger bias (see eq.(9l), 
because this temperature generally overestimates the electron tem- 
perature. 

Including dynamically disturbed clusters does not signifi- 
cantly affect the average value of the recovered Da- However, the 
resulting distribution is clearly asymmetric and presents a large tail 
towards low Da values. In fact, eqs. and J2j show that that 
Da oc (rie)^ / (n^) . Therefore, the presence of clumps in the gas 
distribution produces an underestimate of Da by this factor with 
respect to a completely smooth gas distribution. By looking at the 
distributions of the (3 and rc (central and left panels of Fig|6|, 
unrelaxed structures tend to have rather flat gas density profiles. 
Fitting them with a /3-model forces the slope to be very small, 
with a preference for the core radius to be consistent with zero. 
For instance, the cluster shown in Fig|4| requires /3 = 0.52 and 
rc/ R500 = 0.03, while its estimate of the angular-size distance 
gives DT/DT" = 0.17. 

As a word of caution in interpreting such results, we empha- 
size that this bias in the Da estimate, due to the isothermal gas as- 
sumption, is likely to represent an overestimate of the actual effect 
in real cluster observations for at least two reasons. First, radiative 
simulations of clusters are known to produce temperature gradi- 
ents that, in the central regions, are steeper than observed (see the 
discussion in §2.1). As a consequence, simulated clusters exagger- 
ate the departure from isothermality. Second, the /9-model fitting 
to the Compton-j/ profile has been performed by assigning equal 
weight to all radial bins, with the more external regions bringing 
down the overall normalization of the model profile. In a realistic 
observational setup, the signal from central cluster regions should 
have a relatively larger weight, thus reducing the bias in the recov- 
ered yo. Addressing appropriately this issue would require imple- 
menting detailed mock SZ observations of our simulated clusters, 
a task that we defer to a future analyses. Even keeping in mind 
these warnings, it is clear that deviations from isothermality must 
be taken into account for a precise calibration of the cosmic dis- 
tance scale from the combination of X-ray an d SZ observations of 
galaxy clusters (e.g.. lUdomprasert et all2004l) . 



4.2 Results from the polytropic fit 

In the case of a more general polytropic equation of state, the pa- 
rameters l3 and 7 are calculated by requiring the model to reproduce 
at the same time both the X-ray surface brightness and the temper- 
ature profiles. After obtaining the core radius rc and the normal- 
ization Sxo from the X-ray profile, and To from the temperature 
profile, we combine the two exponents in eqs. and J4j to derive 
both l3 and 7, with yo finally obtained from the SZ profile. 

In FigureQwe show the temperature and Compton-y profiles 
for our example cluster, along with the best-fitting predictions of 
the polytropic /3-model, for the three different definitions of tem- 
perature. The polytropic equation of state provides a reasonable ap- 
proximation to all temperature profiles and, unlike the isothermal 
case, allows us to correctly predict also the Compton-j/ profile. The 
corresponding distributions of /3 and 7 are shown in Figure|8|(we 
do not report the distribution of rc, since it is, by definition, iden- 
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Figure 6. The distribution of the values of /J^sc^^jtrue ^gjj panel), of /3 (central panel) and of rc (in units of -R500; right panel). The dashed and the 
solid lines are respectively for the whole sample and for the subset of clusters classified as regular. Also shown with the light and dark gray areas are the 
con'esponding distributions for the subset of the clusters with 2.5 < T^j (keV) < 5 and T^i (keV) > 5, respectively, classified as regular. The distribution of 
j-)rec ^ £)true jg Q^tained by using the electron-weighted temperatures of the simulated clusters in eg.lgl. 



tical to that of the isothermal model). For both quantities, the ef- 
fect of using different definitions of temperature is rather small. 
As expected, using a polytropic temperature profile implies only a 
modest decrease of the (3 values, because of the weak temperature 
dependence of the cooling function. All the three distributions of 
7 h ave an average value ~ 1.2, similar to observational estimates 
(e.g. lDe Grandi & Molendil2002 h. Moreover, the scatter in this dis- 
tribution is so small to make the isothermal ICM an extremely un- 
likely event. 

The results obtained for Da are shown in Figure |9| and 
also reported in Table Q using emission-weighted, electron and 
spectroscopic-like temperatures. Quite interestingly, the improved 
quality of the fit to the profile of the Compton-y parameter now 
makes the distribution peak at a value much closer to the correct 



Da, independently of whether we use the whole sample or the sub- 
sample of relaxed clusters. 

The angular-diameter distance is correctly recovered when us- 
ing either the electron or the spectroscopic-like temperature with 
deviations which are always i$ 5 per cent, on average. This is a 
rather encouraging result, since it indicates that any bias, induced 
by using the temperature as measured from X-ray observations, is 
in fact rather small. Using instead the emission-weighted temper- 
ature turns into a systematic underestimate of Da by about 20 per 
cent, as a consequence of the fact that it is systematically higher 
than the electron temperature. For all the definitions of tempera- 
ture we find a significant intrinsic scatter, of about 20 per cent on 
average, in spite of our selection of regular objects. 

The fact that the scatter is stable against the definition of tem- 
perature implies that it is almost insensitive to the thermal struc- 
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ture of the ICM and, therefore, to the details of the ICM physics. 
This scatter instead quantifies the effect of departure from spher- 
ical symmetry of the ICM spatial distribution. In fact, the above 
scatter increases to about 50 per cent, if no preselection of regular 
clusters is implemented (Table 0. Quite remarkably, the intrinsic 
scatter calibrated with our simulations is rather close to the 17 per 
cent value, reported by Hughes & Birkinshaw il998h . for the un- 
certainty induced by the intrinsic cluster ellipticity. 

Similarly to the case of the isothermal fit, we note from Fig.|9| 
that the 1ow-Da tails of the distributions are contaminated by irreg- 



ular clusters, for which Da is very badly recovered. For instance, 
for the irreg ular cluster shown in Fig.|4lwe find i3;j7D^™ = 0.30 
when using the electron temperature. Similarly to the case of the 
isothermal fit, also in this case the distribution of the hot clusters is 
consistent with that of the regular cluster subset. 
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4.3 Implications for cosmological parameters 

The precision in the recovery of the angular-size distance when us- 
ing the polytropic model indicates that this method is potentially 
accurate to estimate cosmological parameters. In order to test this 
we create a simple mock catalog of clusters, which is obtained by 
distributing 2/3 of our regular clusters uniformly in redshift in the 
range 0.1 < z < 0.5, while the remaining 1/3 is distributed uni- 



formly in the range 0.5 < z < 1.5. Recall that each simulated clus- 
ter is observed along three orthogonal lines of sight and the redshift 
of each projection is chosen randomly. Figure fTol shows the result- 
ing distribution of clusters in the Da-z plane. We remind here that 
our simulated clusters have been identified at z = 0. Therefore, 
our procedure to distribute them at 2: > neglects the effect of 
their possible morphological evolution. We have been forced to this 
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Table 1. The values of the accuracy in recovering the angular-diameter dis- 
tance, using both the isothermal and the polytropic model, 
and using the emission-weighted, the electron and the spectroscopic-like 
definitions of temperature. For each of them, the first line reports the mean 
and standard deviation, the second the median and the limiting values en- 
compassing 68% of the data. 



choice by the small volume of our simulation box, which implies 
the rapid disappearance of reasonably massive clusters inside the 
high-redshift simulation box. 

For the estimate of the Hubble constant, Hq, we limit the anal- 
ysis to the 66 clusters at 2 < 0.3. Including high redshift ob- 
jects would make the recovery of Hq progressively more depen- 
dent on the knowledge of the underlying cosmology. When using 
the electron temperature, the distribution of the Ho values has mean 
Ho = 70 ± 2 kms~^Mpc^^; when using the spectroscopic tem- 
perature this mean is Ho = 75 ± 2kms~^Mpc~^. In both cases 
the uncertainties are the 1 — cr standard deviations. These values 
are obtained by assuming the correct cosmology. When assuming 
the Einstein-de-Sitter model, we find Ho biased low by 8 per cent. 

As for the estimate of the matter density parameter ilm , we fix 
the value of Ho to its true value and assume flat geometry. In this 
case, we use the 73 clusters lying at z > 0.5. Estimating flm as the 
average of the values yielded by each cluster would provide unre- 
liable results; in fact, inaccurate values of Da can imply negative 
values of ilm , which are clearly unphysical. Therefore, we compute 
the best-fitting value of Qm with a x^-minimization procedure. To 
associate the uncertainty to the estimated Qrn, we reso rt to a boot- 
strap resampling procedure (e.g. §15.6 of IPress et al .' 1992). Each 
bootstrap sample is constructed by randomly selecting, with repeti- 
tion, the objects from the original sample. Each time that a cluster 
is selected, its Da is perturbed with a Gaussian random shift with 
variance 20%, independently of redshift, to account for the "ob- 
servational" uncertainties. The application of this procedure, when 
using the electron temperature, gives $1™, = 0.29±0.05; we obtain 
fim = 0.36±0.06, when using the spectroscopic-like temperature. 
The uncertainties are the 1 — a standard deviations computed with 
100 bootstrap resamplings. The two temperature definitions pro- 
vide two fern's whose difference is consistent with the difference 
in the median values of Da- Moreover, and reassuringly, in both 
cases the central values are consistent with the true value of Qm ■ 

The small size of the errobars of the estimated fim's should 
be clearly taken with caution for at least two reasons. First of all, 
we have assumed errors in Da to be 20 per cent, independently of 
redshift. High-quality SZ and X-ray observations will eventually 
allow to bring statistical errors down to this level in the near future. 
Of course, systematic errors in SZ observations, associated for in- 



stance to point-source contamination and CMB signal removal, are 
different in nature and more difficult to eliminate. 



5 CONCLUSIONS 

In this paper we have applied the method to calibrate the cos- 
mic distance scale from the combination of X-ray and Sunyaev- 
Zeldovich (SZ) observations to an extended set of hydrodynami- 
cal simulations of galaxy clusters. The simulations have been per- 
formed with the GADGET2 code, for a flat ACDM model with 

= 0.3, h — 0.7 and erg = 0.8, and include the effect of 
cooling, star formation and supernova feedback. The aim of our 
analysis was to understand the possible biases introduced by the as- 
sumptions of isothermal gas and the X-ray temperature as a close 
proxy to the electron temperature, as usually done in the analysis of 
real clusters. Furthermore, the application of this method to a large 
set of simulated clusters allows us to quantify the intrinsic scatter 
associated with a cluster-by-cluster variation of their shapes. 

Our main results can be summarized as follows. 

(a) Neglecting the temperature gradients in the application of the 
/3-model produces a significant underestimate of the central value 
of the Comptonization parameter, yo- In turn, this introduces a se- 
vere bias in the estimate of the angular-size distance, Da ■ 

(b) Accounting for the presence of the temperature gradients with 
a polytropic /9-model substantially reduces this bias to a few per 
cent level. While this result holds when using either the electron or 
the spectroscopic-like temperature, using the emission-weighted 
temperature gives a ~ 20 per cent underestimate of Da ■ 

(c) Cluster-by-cluster variations of the asphericity and of the de- 
gree of gas dumpiness cause an intrinsic dispersion of about ~ 20 
per cent in the estimates of Da- This dispersion significantly in- 
creases in case unrelaxed clusters are included in the analysis. 

(d) The set of simulated clusters is used to generate a mock sam- 
ple of clusters out to redshift z — 1.5. By assuming a 20 per 
cent precision in the estimate of Da for each cluster, we find 
that the correct value of Ho is recovered with a statistical error 
of 2kms~^ at la. Furthermore, assuming a prior for the Hubble 
constant and flat geometry, we find that also the matter density pa- 
rameter can be estimated in an unbiased way with a statistical error 
of Af7„ = 0.05. 

It is worth reminding here that our results are based on the 
analysis of simulated X-ray and SZ maps, which are ideal in a num- 
ber of ways. First of all, they have been generated by projecting the 
signal contributed by the gas out to about six virial radii. A more 
rigorous approach would require projecting over the cosmological 
light cone, to properly account for the fore/background contamina- 
tion. While projection effects ought to be marginal f or the X-ray 
maps , they may substant ially affect the SZ signal fe.g.. lWhite et alJ 
l2002HDoIag etall2005h . Furthermore, our noiseless maps need to 
be properly convolved with the "response function" of both X-ray 
and SZ telescopes under realistic observing conditions. Neglect- 
ing the observational noise clearly leads to an underestimate of 
the uncertainties in the determination of the parameters defining 
gas density and temperature profiles. Accounting for such effects 
would definitely require passing our ideal maps through suitable 
tools to sim ulate X-r ay (e.g., Gardini et al. 2004) and SZ obser- 
vations (e.g.. lKneissl e t al. 2001; Pierpaoli et al. 2005). Finally, the 
effect of neglecting the departure from isothermality depends on the 
physical description of the ICM provided by the simulations. Since 
simulated clusters have central temperature gradients, which are 



Angular diameter distance estimates from the SZ effect 1 1 



steeper than the observed ones, the above effect is probably over- 
estimated. This demonstrates that a proper use of hydrodynamical 
simulations to calibrate galaxy clusters as standard rod also requires 
a correct description of the physical properties of the intra-cluster 
gas. 

When this paper was read y for submission it came to our atten- 
tion a paper bv lHallman eio] (2005), which reports on an analysis 
similar to ours. However, their results indicate that assuming an 
isothermal gas tends to overestimate rather underestimate Da, as 
we find. Clearly, a detailed comparison of their approach with ours 
would be necessary to understand this discrepancy. We also note 
that these authors propose a novel method to estimate Da which 
relies on spatially resolved SZ observations. On the contrary, our 
suggestion of considering a polytropic gas model can also be ap- 
plied when the spatial information of the SZ data is poor. 
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